Confinement of electrons in size modulated silicon nanowires 
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Based on first-principles calculations we showed that superlattices of periodically repeated junc- 
tions of hydrogen saturated silicon nanowire segments having different lengths and diameters form 
multiple quantum well structures. The band gap of the superlattice is modulated in real space as 
its diameter does and results in a band gap in momentum space which is different from constituent 
nanowires. Specific electronic states can be confined in either narrow or wide regions of superlattice. 
The type of the band lineup and hence the offsets of valence and conduction bands depend on the 
orientation of the superlattice as well as on the diameters of the constituent segments. Effects of 
the SiH vacancy and substitutional impurities on the electronic and magnetic properties have been 
investigated by carrying out spin-polarized calculations. Substitutional impurities with localized 
states near band edges can make modulation doping possible. Stability of the superlattice structure 
was examined by ab initio molecular dynamics calculations at high temperatures. 

PACS numbers: 73.63.Nm, 73.22.-f, 75.75.+a 



I. INTRODUCTION 

Rod-like Si nanowires (SiNW) have been synthesized 
down to oj 1 nm diameter PJ They are attractive one- 
dimensional (ID) materials because of the well-known 
silicon fabrication technology that make them directly us- 
able on the Si-based chips. Even if unsaturated dangling 
bonds on the outer surface usually attribute a metal- 
lic character to SiNWs, they become insulator (or semi- 
conductor) upon saturation of these dangling bonds by 
hydrogen atoms. 2 SiNWs display diversity of electronic 
properties depending on their diameter, as well as their 
orientation. In particular, the band gap of semiconduc- 
tor SiNWs varies with their diameters. They can be 
used in various electronic, spintronic and optical appli- 
cations, such as field effect transistors^, light emitting 
diodes^, lasers^ and interconnects. The conductance of 
these sem iconductor nanowires can be tuned easily by 
doping 6 7 during the fabrication process or by applying a 
gate voltage. Recent studies have shown that 3d transi- 
tion metal doped Si nanowires become half-metallic. 8 

This letter demonstrates that SiNWs of different diam- 
eters can form stable superlattices. The electronic band 
structure of the superlattice is different from the con- 
stituent SiNWs and is modulated in real space leading to 
a multiple quantum structure and/or to a series of quan- 
tum dots. In these size induced quantum wells, specific 
states are confined. One dimensional multiple quantum 
well structures generated by compositional modulation 
of nanowires were examined previously. For example, 
superlattices of GaAs/GaF^, InAs/InF^ and Si/SiGeP 
nanowires were fabricated. Moreover, superlattices of Si 
and Ge were investigated theoretically.'^ 2 -' In these struc- 
tures multiple quantum wells are formed because of the 
different nature of materials that constitute the nanowire 
superlattice. In the present work however, multiple quan- 
tum wells are formed because of the quantum size effect, 
which is a diversification of the electronic structure of the 



same material with a change in its size. 

Our results are obtained by performing first-principles 
plane wave calculations within Density Functional The- 
ory (DFT)P using ultrasoft pseudopotentialsP^ The 
exchange correlation potential has been approximated 
by Generalized Gradient Approximation using PW91 
functional. A plane-wave basis set with kinetic energy 
cutoff of up to 250 eV has been used. All atomic positions 
and lattice constants are optimized. The convergence for 
energy is chosen as 10~ 5 eV between two steps, and the 
maximum force allowed on each atom is 0.02 eV/A. More 
details about calculations can be found in Rcf[8j. 



II. SIZE MODULATED SILICON NANOWIRES 

Here Si nanowires (i.e. SijVi with diameter ~ D\ 
and having Ni Si atoms in the primitive unit cell) are 
first cut from the ideal bulk crystal along desired direc- 
tion. Then in every alternating segment comprising l\ 
unit cell the diameter D\ is kept fixed, but in the adja- 
cent segment comprising Z 2 unit cell (each having N 2 Si 
atoms) the diameter is reduced to Z?2- The latter part 
can be identified as the segment of Sijv 2 nanowire. At the 
end, the segments of Si and Si n 2 have made a smooth 
junction and hence formed an ideal superlattice Sijv so 
that its diameter is modulated in real space. Note that 
N < hNi + I2N2, because some surface atoms attaching 
with a single bond to the surface were removed at the be- 
ginning. Subsequently we relaxed the atomic structure 
of this bare Si at. Upon relaxation, the dangling bonds 
on the surface are saturated by H atoms and SijvHm 
superlattice is further relaxed for final atomic structure 
and the lattice constant. The resulting superlattice can 
be described by a Si rod with alternating diameters or a 
nanowire with alternating wide and narrow parts. Here 
we consider two superlattice structures, namely Sii3sH42 
and Sii57H64, which are grown in [111] and [100] direc- 
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FIG. 1: Energy band structure of the superlattice, Sii3sH42, 
formed from periodically repeated junctions consisting of one 
unit cell of SiNW with D^llA and two unit cell of SiNW 
with D2 ~ 7 A formed along [111] direction. Here the diame- 
ter D is defined as the largest distance between two Si atoms 
in the same cross-sectional plane. Isosurface charge densi- 
ties of specific bands and their planarly averaged distribution 
along the superlattice axis together with their confinement in 
percentages are also shown. Large and small balls indicate Si 
and H atoms. Zero of energy is set at the Fermi level shown 
by dash-dotted line. 



tions, respectively. For the latter we investigate the sur- 
face defect and also boron, B and phosphorus, P substi- 
tutional impurities. 



A. Sii38H42 Superlattice 

To form the superlattice, Sii3gH42, we took l\ = 1 
unitcell of Si68 and 1? = 2 unitcell of Si^s cut in [111] di- 
rection of the bulk silicon crystal. The combined struc- 
ture had six Si atoms attached to the surface by only 
one bond. These atoms were removed at the beginning 
and the structure was optimized. After the relaxation, 
all Si atoms were at least triply coordinated. The dan- 
gling bonds of triply coordinated Si atoms were satu- 
rated with hydrogen atoms and the structure was relaxed 
again. The band structure and isosurface charge densities 
of the resulting superlattice is shown in Fig. [T| Propa- 
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FIG. 2: Band structures of infinite and periodic Si4sH2o, 
Si25H4 and Si2iHi2 nanowires, which constitute Sii57H64 
structure. The front and side view of these relevant structures 
are shown on the top of the band structure plots. Numerals 
as subscripts of Si and H indicate the number of Si and H 
atoms in the unit cell. Energy band gaps are shaded. 



gating states, with dispersive bands have charge densities 
distributed everywhere in the superlattice rod. For ex- 
ample, a state near ~ -1 eV is a propagating state. The 
states of flat degenerate band at the top of the valence 
band are confined in the narrow segment. The integral 
of the planarly averaged charge, J t \'i'(z)\ 2 dz = 0.96 in- 
dicates a rather strong confinement. The states of the 
lowest conduction band with low dispersion are localized 
in the wide segment with relatively weaker confinement. 
This situation implies the staggered band alignment with 
confined hole states in the narrow part and confined elec- 
trons in the wide part. Confinement may increase with 
the size (length) of the confining segment since the en- 
ergy of the state relative to the barrier varies. For perfect 
localization, i.e. ~ 100 % confinement in l\ or I2, where 
coupling between confined states are hindered, the seg- 
ment l\ or I2 behaves as a quantum dot. Confinement 
and band edge alignment in the present superlattice are 
reminiscent of the 2D pseudomorphic or commensurate 
semiconductor superlatticesP^ 



B. Pristine Sii57H64 Superlattice and Vacancy 

Here we presents the superlattice formation and con- 
fined states in a different structure, namely Sii57H64. To 
form this structure we took li — 2 unitcell of Si45 and 
I2 = 3 unitcell of Si25 cut in [100] direction of the bulk 
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FIG. 3: (a) Electronic structure and isosurface charge densities of selected states for pristine Sii57Hg4 and defected Sii5eH63 
(with one Si-H on the surface removed) nanowire superlattices. Arrows pointing at the same isosurface charge density plot 
indicates that those states are not affected by the formation of the defect and have nearly the same charge density profile. Fermi 
level, represented by dashed-dotted (red) line, was shifted to the valence band edge and set to zero. Solid (blue) and dashed 
(green) lines in the right hand side box represent the majority and minority spin bands, respectively, (b) The distribution of 
interatomic distances of relaxed superlattice. The arrows indicate the ideal bulk 1st, 2nd and 3rd nearest neighbour distances. 



silicon crystal. The combined structure had eight sili- 
con atoms that were attached to the surface by only one 
bond. We have removed these atoms at the beginning 
and relaxed the resulting structure. Then we have pas- 
sivated the dangling bonds of the final structure and re- 
laxed the whole structure again. As a result, the wide 



part of the nanowire superlattice had l\ — 2 unitcell 
of Si45H2o with D\ ~ 14 A and the narrow part had 
h — 3 unitcell of Si25H4 and Si2iHi2 mixture with D2 
~ 7 A. Figure [2] presents the band structure of the cor- 
responding infinite periodic nanowires which constitute 
the Sii57Hg4 superlattice structure. To compare with the 
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FIG. 4: (a) Energy per unitcell of Sii57He4 superlattice as its 
temperature is increased from 500 K to 1400 K (green/light) 
and kept constant at 1300 K (red/dark). Trajectories of atoms 
at the surface of (b) wide and (c) narrow parts and at the 
center of (d) wide and (e) narrow parts of the superlattice, 
projected on a plane perpendicular to the wire axis, as the 
system is kept at 1300 K. 

band structure of the superlattice (presented in Fig. [3] 
as pristine) , the band structure of these nanowires were 
folded l\ + 1% = 5 times. One can see that the band gap 
of the superlattice is close to the lowest band gap of its 
constituents. This is expected for the superlattice struc- 
tures having normal band lineup. The band gaps of struc- 
tures shown in Fig. [2] do not obey the well-known trend, 
Eq = C/D a + E G \bulkSi], occurring due to the quantum 
confinement effectP^ at small D. We ascribe this result 
to the surface reconstructions which differ in accordance 
with the procedure of atomic relaxation. Quantum con- 
finement trend occurs when Si nanowires, cut from ideal 
bulk, are directly passivated with H atoms and subse- 
quently relaxed. In the present study, however, the ideal 
structures are first relaxed, then passivated with H atoms 
and then relaxed again to obtain the final structure!^ 
Such Si nanowire structures with D~l nm may not obey 
the quantum confinement trend. The decline of the gap 
values from the 1/D a behaviour has been also reported 
for other types of SiNWs and was attributed to surface 
effectsP 



Figure |3ja.) presents the band structure and the pro- 
jected charge density isosurface plots of pristine Sii57H64 
and defected Sii56H 63 nanowire superlattices. Flat mini- 
bands of Sii57H64 structure near the edge of conduc- 
tion and valence bands are distinguished. 1 In contrast 
to propagating states of dispersive bands near 1.5 eV, 
the states of these flat mini bands are confined in the 
narrow segments of the superlattice. This situation also 
implies a normal band lineup. We are aware of the fact 
that fabrication of nanowire superlattices having ideally 
tapered as in our case is impossible. So we examined the 
effect of an imperfection on the behavior of the system 
by removal of a Si-H atom pair, which was normally at- 
tached to the surface. The location of this defect and 
structural change after relaxation is given in Figure |3ja) . 
The resulting structure has odd number of electrons, so 
we made spin polarized calculation. One can see that 
the valence and conduction band edges are not affected 
by defect formation, while one empty and one filled state 
localized around the defected area appear near the va- 
lence band edge. The magnetic moment induced by the 
defect formation is 1 (1b- The vacancy states in bulk Si 
are located deep in the band gap. In the superlattice 
structure, this vacancy state is located near the valence 
band, since the defect is formed on the surface at the in- 
terface between two segments having different diameters. 

In Figure |3jb) the distribution of interatomic distances 
shows significant reconstruction and deformation. In par- 
ticular, owing to increased surface/volume ratio the de- 
viation of the second and third nearest neighbour dis- 
tances from the bulk crystal is rather large at the nar- 
row and interface regions. We note that in both short 
periodicity superlattices discussed above the electronic 
structures are significantly different from the constituent 
SiNWs. As k and/or increases, the superlattice effects 
and confined states may become more pronounced, and 
band-offset converges to a well-defined value. When the 
lengths of l\ and h > Xb (Broglie wavelength), the elec- 
tronic properties in each segment becomes close to those 
in constituent nanowires. In addition to l\ and I2, Di 
and D2 and also the direction of growth are also criti- 
cal parameters, since the electronic structure of the con- 
stituent SiNWs is strongly dependent on these parame- 
ters. As Di and D2 increases, the surface effects decrease 
and the density of states becomes more bulk like. The 
discontinuities of valence and conduction bands of the su- 
perlattice is expected to diminish, since the band gap of 
the constituent nanowires converge to the bulk Si value. 
However for a significant AD = D\ — D2 the superlattice 
effects, in particular confinement of specific states are ex- 
pected to continue. In contrast, smooth interfaces lead- 
ing to hornlike profiles can be used for focused electron 
emission.^ Recently, transport properties of size modu- 
lated SiNWs have been shown to depend on the growth 
directionPS 

The stability of the Sii57H64 structure was analyzed us- 
ing the finite temperature ab initio molecular dynamics 
calculations with time steps of 2 x 10~ 15 seconds. Fig- 




FIG. 5: Atomic and electronic structure and charge density isosurfaces of Sii57H64 substitutional doped by B and P atoms 
(namely Sii56P(B)H64). Top panel presents a view of the superlattice structure, where the doped site is represented by 
large red(dark) ball. Arrows between the bands of the undoped (pristine) superlattice and the bands of the B and P doped 
superlattices indicate the states, which preserved their character even after doping. The charge density isosurfaces of specific 
states are linked to the corresponding bands by arrows. Fermi level, represented by dashed-dotted (red) line, was shifted to 
the valence band edge and set to zero. Solid (blue) and dashed (green) lines represent the majority and minority spin states, 
respectively. Small panel at the left hand side presents the charge density isosurface plot of B-Si bonds. 



ure 4ja) presents the total energy difference of the system 
with respect to the ground state energy at K as a func- 
tion of time. In the first calculation we have raised the 
temperature of the system from 500 K to 1400 K in 2 ps. 
During this calculation there was no major change in the 
structure and the interface region was not destroyed by 
diffusion. In the second calculation we took atomic posi- 
tions and velocities of the system from the first calcula- 
tion when the system was at 1300 K and kept the tem- 
perature constant for 9 ps. This calculation also resulted 
in no major change in the structure of the superlattice. 



Figure Qb) and (c) present the path of Si atoms during 
the second calculation at the surfaces of wide and narrow 
parts respectively. The range of the thermal fluctuations 
are nearly the same and span a diameter of about 1 A. 
Figure Qd) and (e) present the path of Si atoms at the 
center of the wide and narrow parts respectively. As 
seen the range of thermal fluctuations in the narrow part 
(~ 0.9A) is larger than that in the wide part (~ 0.5A). 
It should be noted that, the temporal trajectories of 9 
picoseconds in the present ab initio molecular dynamics 
calculations are rather long as compared to ones usually 
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carried out in the stability tests, but may not be long 
enough to accumulate necessary statistics. However, the 
temporal trajectories at temperature as high as 1300 K 
in the present study are long enough to eliminate the ex- 
istence of any atomic configuration in a weak local min- 
imum of the Born-Oppenheimer surface, which is prone 
to structural instability. Therefore, based on the accu- 
rate structure optimization and temporal trajectories as 
long as 9 picoseconds carried out at 1300 K, we believe 
that superlattice structure under study is stable at least 
at a temperature higher than the usual operation tem- 
perature (~425 K) of a device. 



C. B and P doping of Sii 5 7He4 Superlattice 

Next we examine the effect of doping of the Sii57Hg4 
superlattice by substitutional B and P impurities. Both 
structures have odd number of electrons and hence spin 
polarized calculations have been carried out. Neverthe- 
less, B doped structure resulted in a nonmagnetic state 
with half filled band edge. P doped structure resulted 
in a magnetic state with a magnetic moment of 1/i_b. In 
Fig. [5] the electronic structure of pristine and doped su- 
perlattices are shown to follow the effect of the doping 
on the states at the band edges. Owing to the small di- 
ameter of the superlattice, the band structure is affected 
by the impurity. Since the band gap of the wide part 
is larger than that of the narrow part, the state asso- 
ciated with the impurity, which occurs normally either 
near the valence or conduction band, falls in the bands 
of the narrow part. 

It is known that substitutional B doping of bulk silicon 
crystal results in an acceptor state ~45 meV above the 
valence band edgeP^l In the substitutional doping of the 
superlattice by B similar situation occurs and the accep- 
tor state is located above the valence band edge corre- 
sponding to the wide part, but falls in the valence band 
of the narrow part. The filled first and the lower lying 
three valence bands at about ~-0.4 eV undergo changes: 
The filled first valance band having states confined in 
the narrow part becomes half-filled. Three lower lying 
valence bands of the pristine superlattice are raised to- 
wards the edge of the valence band edge. The states 
of these bands mix with the tetrahedrally coordinated 
sp 3 hybridized orbitals of B and form a triply degener- 
ate state, just below the half filled valence band edge. 
Charge density isosurfaces of these states around boron 
atom presented in Fig. [5] clearly indicates sp 3 hybridiza- 
tion in tetrahedral directions. As a result, the "hole like 
states" becomes confined in the narrow part. As for the 
conduction band is slightly shifted down, but its states 
preserve their character. 

In P doping of the superlattice at the center of the 
wide part, the impurity state localized at P atom occurs 
in the conduction band of the superlattice. This state, 
in fact, the donor state occurring ~ 200 meV below the 
conduction band corresponding to the wide part. This 



energy is in good agreement with the energy calculated 
recently for the ionization energy of P donor states in 
the Si nanowires.^ Below this donor state, all conduc- 
tion band states are confined cither in the narrow part or 
at the interface, but not in the wide part. Energy shifts 
of the band edge states are more pronounced in structure 
having substitutional P doping. Conduction band edge is 
shifted down by 0.75 eV, but the charge density profile re- 
mains the same. One spin-up and one spin-down bands 
split by ~0.2 eV are formed near the conduction band 
edge with charge densities localized in the narrow part. 
The triply degenerate state in the valance band near the 
edge are shifted down. That states localized around sub- 
stitutional P impurity appear in the conduction band 
continua, while a filled spin up band states are localized 
in the defect free narrow part, indicate the possibility of 
modulation doping. In closing this section we note that 
the electronic structure associated with the doping of the 
small diameter superlattices is complex and has aspects 
different from that in the doping of bulk Si. The small 
size and dimensionality of the doped superlattices hin- 
ders the application of the effective mass approximation 
leading to hydrogenic impurity model. For the same rea- 
son, the relative positions of band edges can be affected 
by the presense of a substitutional impurity. 



III. CONCLUSION 

In conclusion, we considered two specific superlattices 
grown in different crystallographic directions as a proof of 
concept. Each superlattice can be viewed to form the pe- 
riodically repeating junctions of two hydrogen saturated 
silicon nanowires of different diameters. Since the band 
gaps of the constituent Si nanowires are different, the 
band gaps of these superlattices are modulated in real 
space. As a result, one dimensional multiple quantum 
well structures with states confined in different (wide or 
narrow) parts are formed. At strong confinement or lo- 
calization, the superlattice structures can be viewed as a 
series of quantum dots. We believe that the bad gap 
modulation and resulting quantum well structure and 
confinement can occur as long as the band gap differences 
are maintained between adjacent region at both sides of 
the junction. Furthermore, we showed that the electronic 
structure resulting from the substitutional doping of the 
superlattice interesting and give rise to modulation dop- 
ing. We finally note that, the band gaps of SiNWs at 
both sides of junction forming the superlattice structure 
are underestimated by DFT/GGA calculations. Since 
the systems under study are too large, these bands can- 
not be improved by GW calculations. Nevertheless, our 
work is just a proof of concept. According to our results, 
if the band gaps are different in different parts of the het- 
erostructure, confined states can occur and modulation 
doping of these quantum structures can be achieved. De- 
tails of energy band structure and location of impurity 
states in the gap should only be taken qualitatively. 
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Various types of electronic devices, such as resonant 
tunnelling double barriers generated from these super- 
lattices, can be arranged on a single rod, where they can 
be connected by the metallic segments of unsaturated Si 
nanowires. The character of these devices can be engi- 
neered by varying their growth direction and structural 
parameters, such as 1%, l 2l D\, D 2 - Even if the diameters 
of superlattices we treated in this study are small (~ 1 
nm) and hence in the range where band gaps vary signif- 
icantly with diameter, they are still in the range of those 



fabricated recently.^ 
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